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Abstract 

An asymmetric generalization of the zero-temperature g-state Potts 
model on a one dimensional lattice, with and without boundaries, has 
been studied. The dynamics of the particle number, and specially the large 
time behavior of the system has been analyzed. In the thermodynamic 
limit, the system exhibits two kinds of phase transitions, a static and a 
dynamic phase transition. 



1 Introduction 



In recent years, reaction-diffusion systems have been studied by many people. As 
mean-field techniques, generally, do not give correct results for low-dimensional 
systems, people are motivated to study exactly-solvable stochastic models in 
low dimensions. Moreover, solving one-dimensional systems should in princi- 
ple be easier. Exact results for some models on a one-dimensional lattice have 
been obtained, for example in [p]-|l2|, and in Jll|,[l4| for the case of multi-species 
systems. Different methods have been used to study these models, including an- 
alytical and asymptotic methods, mean field methods, and large-scale numerical 
methods. 

Some interesting problems in non-equilibrium systems are non-equilibrium 
phase transitions described by phenomenological rate equations, and the way the 
system relaxes to its steady state. Kinetic generalizations of the Ising model, for 
example the Glauber model or the Kawasaki model, are such phenomenological 
models and have been studied extensively |l5|~|l1| . Combination of the Glauber 
and the Kawasaki dynamics has been also considered [pO| [2^] . One of the models 
which has been extensively studied is the g-state Potts model evolving according 
to some generalization of the zero-temperature Glauber dynamics (see p3|-p8| , 
for example). The evolution tends to align all the spins, and so domains of 
parallel spins grow with time. 

In this paper, we want to study an asymmetric generalization of the zero- 
temperature q-state Potts model on an infinite lattice with and without bound- 
aries. There are also reactions at the boundaries. The dynamics of the particle 
number, and specially the large time behavior of the system is studied. In the 
thermodynamic limit, the system shows two kinds of phase transitions. One of 
these is a static phase transition, the other a dynamic one. The static phase 
transition is controlled by the reaction rates, and is a discontinuous change 
of the behavior of the derivative of the stationary particle density at the end 
points, with respect to the reaction rates. The dynamic phase transition is con- 
trolled by the reaction rates at boundaries, and is a discontinuous change of 
the relaxation time towards the stationary configuration. There is a fast phase 
and a slow phase. Increasing the reaction rates at the boundaries, the system 
goes from the slow phase to the fast phase. This model may be considered as a 
biased voting model, in the sense that there are q different opinions, and when 
people at the adjacent sites have different opinions, they may interact so that 
their opinions become the same. 



2 Asymmetric g-state Potts model at zero tem- 
perature 

In the ordinary Glauber model, the interaction is between three neighboring 
sites. Spin flip brings the system to equilibrium with a heat bath at temperature 
T. A spin is flipped with the rate \x := 1 — tanh[J/(fcT)], if the spin of both 
of its neighboring sites are the same as itself; and is flipped with the rate A := 



f 



1 + tanh[J/(fcT)], if the spin of both of its neighboring sites are opposite to it. 
At domain boundaries, the spins are flipped with unit rate. So the interactions 
can be written as, 

AAA — > A$A and 000 -> 0A0 fx 

A§A — > AAA and 0/10 -> 000 A 

AA0^A00 and 00A ^ 0/L4 1 

where spin up and spin down are denoted by A and 0. One can interpret an up 
spin as a particle, and a down spin as a hole. At zero temperature, the Glauber 
dynamics is effectively a two-site interaction |27]] : 

A% -> (AA, 00) U (AA, 00) (1) 

where all the above processes occur with the same rate. 

One can consider the following interactions, as an asymmetric generalization 
of the zero-tempcrature Glauber model. 



AA, 



and 



AA 



v 



U^{Z (3) 



u 



If u 7^ v, the above system has left- right asymmetry. The above system on 
an infinite lattice has been investigated in [^9| , where its n-point functions, its 
equilibrium states, and its relaxation towards these states are studied. It can be 
easily shown that the time evolution equation for the average particle-numbers 
of the system with the above interactions are the same as that of a system with 
the following interactions, where diffusion is also present: 



(4) 





'04 


A 


A0^ < 


AA, 


u 




00, 


V 



and 



%A 



A0, n 

AA, v - n (5) 
00, u-fi 



The model addressed in this article is an asymmetric generalization of the 
zero-temperature g-state Potts model in one dimension and with Glauber dy- 
namics. The interaction is a nearest neighbor interaction and is defined as 
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follows 

^^\ A ; A r u (6) 

[ApAp, v 

Similar to (0) and (||) , adding diffusion do not alter the time evolution equation 
of the average particle numbers. The symmetric case (it = v) has been recently 
analyzed using a generalization of the empty interval method p^ j. The one- 
species case (q — 2) has been studied in |29| (without boundaries) and Q 
(with boundaries). It is shown that for the infinite lattice without boundaries, 
the system has two ground state; The sites are all occupied or all vacant. This 
is reminiscent of the ordinary Glauber model, and asymmetry does not change 
the ground states. Denoting by |0) (|^)) the empty (full) lattice, the state of 
the system at infinitely large times is 

|P(oo)) = (l-po)|0>+po|n). (7) 

where po is the initial average density. For the asymmetric generalization of the 
g-state Potts model, there are q ground states. To see this, one can divide the 
state at each site to two classes, A\ and not-^i, for example. The state not-Ai 
is denoted by B\. So the interactions are 




and 




Similar to the asymmetric zero-temperature Glauber model, the ground states 
are the states that all the sites of the lattice are in the state A\, or none of the 
sites are in the state A\. Repeating the same argument for the substates of B\, 
it is seen that there are q ground states. Each ground state is a state in which 
the states of all lattice sites are the same state, A a for example. So the final 
state is 

|P(oc))=5>«LA Q ), (10) 

a 

where \ A a ) is the state for which all the sites are in the state A a , and pg is the 
initial average density of the state A a . Then, all the correlation functions at 
infinitely large times can be obtained easily: 

(nfn^--nl) = 6^6^---p^ (11) 

The evolution equation for the average number of the state a at the site j, 

(nf) is 

(nf) = -(« + v)(nf) + u<^_ x ) + v(nf +1 ) (12) 
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The above equation shows that the evolution equations for the average number of 
different states are decoupled, and so the average number of each state depends 
only on the initial value of the average number of that state. In fact, 

(nf(t)) = e-^Y. Im-,(2V^t)<n?(0)>, (13) 

m 

and at large times, 

/v \ i P [-(«+")+2vVl« 

KW)-^0 7t ■ (1 4 ) 

It is seen from the above equation that if u < v, the expectation at the rightmost 
sites tend rapidly to their final value, and obviously for u > v the leftmost sites 
arrive earlier to their final states. 

3 g-state Potts model on a lattice with the bound- 
aries 

In this section, we study q-state Potts model on a lattice with reaction at the 
boundaries. The interactions on the bulk of lattice are (|J). The exchange of the 
states at the first site is 

Ap -> A a with the rate Ag, (15) 

and at the final site, it is 

Ap -> A a with the rate r|. (16) 

For a (3, Ajg and are rates and should be nonnegative. The diagonal 
elements of A and T are chosen so that 

s a k a p = s a Y a p = 0, (17) 

where 

s a = 1. (18) 

From now on, a repeated subscript and superscript imply a summation over the 
repeated index. 

The equations of motion for the average numbers are 



(n?) = -(« + v)(nf) + uinf^) + v(nf +1 ), j^l,L 
(nf) = A|<nf)-«(n?)+«<n?>, 

(h<£) = r%{nl}-u{nZ)+u{nt_ 1 ). (19) 
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Defining the vector Nk through 

^ Q := «> 

equation (|n]) takes the form 

Nk = -(M + v)N k + uN k -i+vN k+1 , k^l,L 
Ni = ATVi - vNi + vN 2 , 
Nl = TN L - uN L + uN L - X . 

For the stationary state, 

Nk = 0. 

Substituting the ansatz 

N k = Bz\ + Cz k 2 
in the first equation of (^l|), one arrives at 

— {u + v) + vzi H = 0, 



(20) 



(21) 
(22) 
(23) 
(24) 



whose solutions are Z\ = 1, z 2 = u/v. First assume u < v. Then, one can 
write (p3|) as 



N k = B + C (~) 



fc-i 



Since s a N^ = 1, it is seen that, B, and C satisfy 

s a B a = 1, 
s a C' a = 0. 



Using the second and third equations of (Plf), one arrives at 



(r-u) 



(A -v)[B + C] + v 

L-l 

J 



s + c'(- 

1> 



B + C'l- 



B + C - 



L-2 



= o, 

= 0. 



In the thermodynamic limit [L — > oo), these two equations yield 

= 0, 

(A-u + n)C" = -AS. 



(25) 



(26) 



(27) 



(28) 



The first equation in ( |28| ) can be used to determine £?. Substituting it in the 
second one, C is obtained. The degeneracy of the zero eigenvalue of T is equal 
to the degeneracy of the stationary average particle number. The matrix T does 
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have a zero eigenvalue, since s is the left eigenvector of T with the eigenvalue 
zero. If the zero eigenvalue of T is degenerate, then the final value of the average 
particle numbers depend on the initial conditions. Also note that the second 
equation of (^8|) has one and only one solution for C . The reason is that the 
real part of the eigenvalues of A are nonpositive. So, v — u is not an eigenvalue 
of A, and hence the matrix A — v + u is nonsingular. 

If u > v, then the term (u/v) k diverges for the right end sites. So, we choose 



„,.... k-L 

ill 



N k = B + C" (J) . (29) 

In the thermodynamics limit, one arrives at 

KB = 0, 

(r + v - u)C" = -TB. (30) 

If u > v, it is the degeneracy of the zero eigenvalue of A, which determines 
the degeneracy of the stationary states of the particle numbers. Here too, the 
second equation of ( |30"1) determines C" uniquely. 

If u < v, the profile of average particle number is flat for the rightmost 
sites. As u exceeds v, it acquires a finite slope, proportional to ln(u/i>). For the 
leftmost sites, this behavior is reversed. In fact, there is a discontinuous change 
of the behavior of the derivative of the stationary particle number at the end 
points. This is a static phase transition, and it is controlled only by the reaction 
rates and does not depend on the reaction rates at the boundaries. The case 
q = 2 has been already studied in 130] . One notes that the matrices A and 7 are 



A = 



—a a 
a' —a 



in terms of the appropriate parameters of that references. Here, the zero eigen- 
values of A and T are nondegenerate (unless the rates are zero). 

Now let's consider the relaxation of the system towards its stationary state. 
The homogeneous part of ( |l9| ) is 

SN£ = h%5N?, (32) 

where 8N is the difference of N from the its stationary value, so that 

s a 5N£ = 0. (33) 

To study the relaxation of the system, we investigate the eigenvalues of h. One 
finds 

Ex%=-{u + v)xt+uxt_ x +vxt +x , k^l,L 
Exl =A%xf -vx% +vx%, 

Ext=T a p x p L -uxt + uxt_ x , (34) 



6 



where the eigenvalue and eigenvector have been denoted by E and x, respec- 
tively. The solution to these equations is 

x%=B a z , t + C a z%, (35) 

where z^s satisfy 

u 

E = -{u + v)+vzi + — . (36) 
z% 

Performing the change of variable Zi —: y^u/vZi, the second and third equations 
of (|34j) take the form 

^E^BZ + CZ- 1 ) ={A-v) s J^{BZ + CZ~ l ) + v (BZ 2 +CZ- 2 ) 

v E(BZ L + CZ- L ) =(T-u)^^j {BZ L + CZ- L ) 

L-l 



(BZ + CZ j (37) 

Substituting E = —(u + v) + ^Jwv(Z + Z^ 1 ) in the above equations, one can 
write them in the matrix form 



-(u + A)Z + -(u + K)Z- 1 + ^/wv \ /s^ 

-(v + r)z L + v ^z L+1 -(v + r)z- L + ^z-l- 1 ) \c 



= 0. (38) 



To have nontrivial solutions for B and C, the determinant of the matrix should 
be zero: 

-(u + A)Z + ^uv -(u + A)Z~ 1 + y^UV \ 

det | = 0. (39) 

-(v + t)z l + ^idz L+1 -(v + r)z- L + ^z- L - 1 ) 



In the thermodynamic limit (L — > oo), if all of the roots of ( p9[ ) are phases, the 
largest real part of the eigenvalues of the Hamiltonian will be — (u + v) + 
from which the relaxation time of the system is 

T f = (u + v - 2 y /mjy 1 . (40) 

We call this the fast phase. 

If some of the roots of (|io| ) are not phases, however, this may be not the 
case: there may be an eigenvalue of the Hamiltonian with a larger real part, 
which corresponds to a larger relaxation time. For Z = r e lS , one has 



Sft(£!) = -(it + v) + y/uv + -^j cost 



(41) 
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This leads to a larger relaxation time, provided 



1 

r + - 

r 



cos(9>2. (42) 



If there is such a solution for Z, the system is said to be in the slow phase. 

So, if none of the roots of (|39| ) satisfy (|4^), then the system is in the fast 
phase, and the relaxation time does not depend on the reaction matrices T and 
A. Otherwise, the relaxation time does depend on the reaction matrices (the 
slow phase). The transition between these two phases is the dynamical phase 
transition. 

Now, let's seek the nonphase solutions of (^9|). If Z is a solution to (gg), l/Z 
is another solution to it. So it is sufficient to seek the solutions with \Z\ > 1. 
In the thermodynamic limit, and for \Z\ > 1, ( |39| ) is simplified to 

det[y/uvZ -(v + T)] det[y^Z -(u + A)] = 0, (43) 

which consists of the characteristic equations for T and A. Denoting the eigen- 
values of F and A by 7 and A, respectively, one has 

Z= V 42, or Z=H±^. (44) 



As r and A are stochastic matrices, their eigenvalues have nonpositive real- 
parts. For u > v, the real part of (v + j)/y/uv is then not greater than 1. Hence 
it cannot satisfy (|42]). So, the only relevant equation for finding the system in 
the slow phase is 

- U + X (45) 

A similar argument shows that for v > u, the first equation of ( |44|) is relevant. 
One also concludes that for u — v, the system has no slow phase. 

So, without loss of generality, let's take u < v. As T is a stochastic matrix, 
it has at least one zero eigenvalue. However, if this eigenvalue is nondegenerate, 
the right eigenvector of T (B) cannot satisfy 

s a B a = 0, (46) 

since s is the left eigenvector of T corresponding to the same eigenvalue. But 
from (|3^), it is seen that B should satisfy (f46|). So, from the eigenvalues of T, one 
should set aside one zero eigenvalue, and consider only the other eigenvalues. 

The system undergoes a dynamic phase-transition at the point that for one 
of the Z's in ( |45| ) the criterion ( f42| ) is satisfied. At this point, the real- and 
imaginary-parts of Z satisfy 



V = ±(.V- 11^/^, X>1, (47) 
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where 



X : 



Y : 



9f(Z). 



(48) 



One can translate this in terms of A. The criterion for the slow phase is then 
seen to be 



A simple way to induce the phase transition is to multiply the matrix A by 
a parameter r. This means multiplying the rates of the reaction at the first site 
by r. As 3? (A) < 0, one can see that for a large enough value of r, the value of 
3?(A) + u — y/wu will be negative (provided 3?(A) ^ 0, that is, provided the zero 
eigenvalue of the matrix A is not degenerate) . So the system will be in the fast 
phase. It is also seen that as r tends to zero, either 2y/wv — -ft(A) — u becomes 
negative, or in the first inequality in (^) the right-hand becomes greater than 
the left-hand side (which tends to zero). So, the system will be in the slow 
phase. Roughly speaking, increasing the reaction rates brings the system from 
the slow phase to the fast phase. A similar argument holds, of course, for the 
case v > u and the eigenvalues of the matrix T. 
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